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ABSTRACT 

The nature of the excess near-infrared emission associated with the magnetic 
white dwarf commonly known as SDSS 1212 is investigated primarily through 
spectroscopy, and also via photometry. The inferred low mass secondary in this 
system has been previously detected by the emission and variation of Ha, and the 
1—2.5 jttm spectral data presented here are consistent with the presence of a late 
L or early T dwarf. The excess flux seen beyond 1.5 /iin in the phase-averaged 
spectrum is adequately modeled with an L8 dwarf substellar companion and cy- 
clotron emission in a 7 MG magnetic field. This interesting system manifests 
several observational properties typical of polars, and is most likely an old inter- 
acting binary with a magnetic white dwarf and a substellar donor in an extended 
low state. 

Subject headings: binaries: close — infrared: stars — stars: fundamental parameters — 
stars: individual (SDSS J121209. 31+013627. 7) — stars: low-mass, brown dwarfs- 
novae, cataclysmic variables — stars: evolution — stars: formation — white dwarfs 

1. INTRODUCTION 

Cataclysmic variables are spectrally exotic binaries, displaying a diversity and com- 
plexity of features, behavior, and variability, effectively masking the spectral signatures of 

1 Gemini Observatory, Northern Operations, 670 North A'ohoku Place, Hilo, HI 96720 

2 Department of Physics & Astronomy, University of California, 430 Portola Plaza, Los Angeles, CA 
90095; jfarihi@astro.ucla.edu 

3 Department of Physics & Astronomy, University of Leicester, Leicester LEI 7RH, UK; mbu@star.le.ac.uk, 

4 Spitzer Science Center, California Institute of Technology, MS 220-6, Pasadena, CA 91125; 
hoard@ipac.caltech.edu 



-2- 



their coolest secondaries. Despite unyielding efforts to unambiguously and directly detect 
likely substellar companions in either non-magneti c or magnetic cataclysmic variables , these 
altered degenerate dwarfs remain largely elusive (ILittlefair. Dhillon. fc Martini 120031 ). The 
ma ss determinatio n of th e subs tellar donors in t wo no n-magnetic accreting binaries reported 
by ILittlefair et al.l (120061 ) and ILittlefair et al.l (120071 ) remain rare among cataclysmic vari- 
ables, although the determination of their substellarity was indirect. Cataclysmic variables 
containing strongly magnetic white dwarfs (called polars because of the highly polarized 
nature of their light) offer a slight advantage by lacking the often overwhelming luminosity 
from an accretion disk. On the other hand, their magnetic fields generate copious cyclotron 
radiation (in addition to flux from direct accretion), which is not present in non- magnetic 
cataclysmic variables. 

Polars are also known for exhibiting long periods of quiescence in whic h many secon- 
darie s have been directly observed, although none of spectral type later than M ([Harrison et al 
20051 ). The recently recognized class of low accretion rate polars may provide decent infrared 
hunting grounds for substellar companions, if present, as their cyclotron harmonics appear 
mainly in the optical due to B ~ 50 MG fields. It has been suggested that no mass transfer 
takes place in these systems, but rather the efficient capture o f secondary wind, an d that 
their separations are consistent with being essentially detached (jSchmidt et al.ll2005bl ). Yet, 
it is difficult to certi fy that these binarie s are not simply in p rotracted low states such as 
seen in other polars (jSchwope et al.l 12007c lHarrison et al.l 120041 ) . 



Recent Spitzer IRAC observations of several polars indicate flat infrared spectral en- 
ergy distributions which are difficult to reconcile with cool degenerate companions alone, 
and may indicate circumbinary dust as w ell as underlying cyclotron continuum emission 
( iBrinkworth et al.ll2007l ; iHowell et al.ll2006l ). It will be interesting to see the results of simi- 
lar mid-infrared observations with a larger sample of polars as well as observations at longer 
wavelengths where there should be little or no cyclotron emission even in B w 10 MG 
fields. The first published Spitzer IRS spectrum of a polar strongly supports the presence of 
circumbinary dust at EF Eridani (jrloard et al. 2007 ). 



This paper represents a spectroscopic effort to directly detect the almost certain substel- 
lar companion to the cool magnetic white dwarf SDSS J121209.31+013627.7 (SDSS 1212). 
Its binarity was reported with optical spectra supplemented with a single J-band photo- 
metric observation, limiting the spectral type of the secondary to L5 or later, and was 
initially thought to represent a detached system (jSchmidt et al.ll2005al ). This possibility was 
still considered viable despite a second study which detected strong cyc lotron emiss i on - a 
likely indicator of mas s tran sfer - in if-band tim e series photometry (jDebes et al. 2006 ). 
Both iKoen Maxtedl (120061 ) and iBurleigh et al.l (j2006al ) detected optical variability and 
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a light curve similar to that shown by polars in low states. iBurleigh et al.l (l2006al ) further 
showed that the amplitude of the modulations was larger at blue and ultraviolet wavelengths 
than in the red , consis tent with a hot accretion spot responsible for the cyclotron emission. 



Burleigh et al.l (j2006al ) also detected x-ray emission from SDSS 1212, conclusively proving 
accretion is taking place. This paper presents cross-dispersed 1 — 2.5 fim spectroscopy and 
spectral modeling, along with HK photometry, in order to constrain the source (s) of the 
observed near-infrared excess emission. 



2. OBSERVATIONS 



2.1. HK Photometry 



Near-infrared images of SDSS 1212 were acquired on 2006 June 10 at the United King- 
dom Infrared Telescope under photometric con ditions (secz < 1 . 1) with the Wide Field 
Camera (WFCAM; Casali et al., in preparation; iHirst et al.ll2006l ; iHenry et al.ll2003l ). A 5 
point dither pattern was used with 2 x 10 s exposures at each position; this pattern was 
repeated twice at H and thrice at K for total exposure times of 200 and 300 s, respectively. 
The data were reduced in the standard fashion. Each stack was median combined to ex- 
tract a sky frame which was also normalized to generate a flat field map. Individual sky 
subtracted, flat fielded frames were registered and averaged for a final frame upon which to 
execute measurements. The full widths at half maximum of the reduced stellar images were 
2.0 and 1.9 WFCAM pixels (0".40 pixel -1 ) at H and K respectively. Aperture photometry 
was performed with standard IRAF tasks using an r = 1".20 aperture and an r = 8 — 12" 
sky annulus. 

Due to the very wide field of view of WFCAM (819" x 819"), 15 2MASS point sources 
in the same field as SDSS 1212 were simultaneously imaged to choose amongst as possible 
photometric calibrators. Each 2MASS star chosen as a calibrator had all of the following 
properties: a full width at half maximum consistent with linear accumulation (i.e. unsatu- 
rated), high signal-to- noise ratios (S/N > 100) in the reduced WFCA M images, and reliabl e 
2MASS photometry at both wavelengths (an < 3.4%, a Ks < 6.3%; ISkrutskie et al.ll2006l ). 
The selected 12 calibrators are listed in Table (U along with their 2MASS data and their 
derived instrumental zeropoints. 

Additional photometry was acquired on 3 July 2006 at the Ge mini Telescope Nort h 
under photometric conditions with the Near-Infrared Imager (NIRI 



rlodapp et all 120031 ). 



Both SDSS 1212 and the UKIRT standard star FS 132 (lHawarden et al.ll200ll ) were observed 
in HK at nearly identical average airmasses of sec z ~ 1.95, yielding extinction corrections 
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less than 0.003 mag. A 5 point dither pattern was used for both science target and calibrator, 
with an extra frame taken for SDSS 1212. The calibrator was observed 10 x 1 s at each 
position while the science target was observed 3 x 12 s at each position, yielding total 
exposure times of 50 s for FS 132 and 216 s for SDSS 1212 at HK. The data were reduced 
in a manner similar to the WFCAM data, with separately acquired and processed flat field 
images. The full widths at half maximum of the reduced stellar images were between 4.5 
and 4.9 NIRI pixels (0".12 pixel -1 at f/6). Aperture photometry was performed similarly 
using an r = 0".47 aperture and an r = 2". 3 — 3". 5 sky annulus. 

The photometry results are listed in Table [2] and plotted in Figure [TJ The third column 
in the table gives the uncertainty in the measured aperture photometry (due to intrinsic 
S/N limitations). The fourth column in the table gives the calibration uncertainty: for the 
WFCAM data, this is the standard deviation in the derived zeropoints; for the NIRI data, 
an intrinsic 3% error was assumed. 



2.2. 1 — 2.5 /im Spectroscopy 

Spectroscopy was performed on 2006 June 19 and 23 at the Gemini Telescope South with 



the Near- Infrared Spectrograph (GNIRS; lElias et al.lll998l ). Observations were executed in 



cross-dispersed mode with the 31.7 1/mm grating and the 0".3 wide slit, yielding R fa 1700 
while effectively covering the entire near-infrared range without inter-order contamination. 
Data were taken at two nod positions along the 6" long slit, with 300 s exposures at each 
nod position, repeated 6 times for a total on source integration time of 1 hr on each of two 
nights. The telluric standard HIP 60030 (A7V, V — K = 0.5) was observed in a similar 
manner with 6 s exposures at each nod position, repeated 4 times. 

The raw data were processed with standard GNIRS tools within the Gemini IRAF 
package version 1.9. These steps included flat field and pinhole map creation from calibra- 
tion lamps, bad pixel correction, linearity correction, sky subtraction, image shifting and 
averaging, order cutting, spatial rectification, and wavelength calibration from argon lamps. 
Spectral extraction of each order was performed with a GNIRS equivalent to the IRAF a pa 1 1 
task. 

Some processing of the raw data was required prior to application of the standard 
package tools listed above. The first and second science frames on each night contained 
significant residual signal (between 3 and 5 times larger than the read level in the unexposed 
portion of the array) from the target acquisition mirror. The raw acquisition images were 
scaled to match the contamination level of the affected frames and subtracted off, effectively 
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solving this problem. Additionally, there is an occasionally occurring bias level problem in 
one quadrant of the detector which appeared in many of the science exposures. This noise 
was well removed with the nvnoise task within the GNIRS tools. 

Both the telluric standard and SDSS 1212 were processed in an identical fashion until the 
spectral extraction step. Instead of the 8 — 10 pixel diameter (1".20 — 1".50) aperture used for 
HIP 60030, a 4—5 pixel diameter (0".60— 0".75) aperture was used for SDSS 1212 to maximize 
the extracted spectroscopic S/N of this faint source. Flux calibration and telluric correction 



were performed for each order w ith Spextool (version 3.3; ICushing. Vacca. fc Raynerl 12004 



Vacca. Cushing. &: Rayner 20031 ). using the general xtellcor package. The flux levels of all 



four spectral orders presented here agree quite well at their boundaries and hence no inter- 
order adjustment was necessary. However, the entire cross-dispersed spectrum, the average 
of both nights' data, was adjusted slightly by a single scaling factor to best match the J 
and //-band photometry in Figure [TJ This unsmoothed spectrum is displayed in Figure [2] 
together with the JHK points from Figure [TJ 



RESULTS 



When telescope time was requested for this study, it was not yet known that SDSS 1212 
was variabl e in any broad wavele ngth bands. Apart from the Ho; line strength variability 
reported by lSchmidt et al.l (j2005al ). SDSS 1212 is no w known to be variable in uqi, VR (and 



unfiltered white optical light), and K, but not in JH (IBurleigh et al.ll2006al ; iKoen fc Maxted 



?r 

20061 : iDebes et alJl200fl ). 



Figure [3] shows the corresponding phase of SDSS 1212 for all the photometric and 
spectroscopic observations. The WFCAM photometry was taken during K-band minimum, 
by which it is meant that all the 2.2 /xm flux should be near its minimum value - whether 
this minimum flux comes from cyclotron continuum emission, the substellar companion or 
a combination of sources is uncertain at present and the focus of this study. The NIRI 
photometry, on the other hand, was taken close to the K-b&nd maximum which occurs near 
orbital phases <p = 0.8 and 0.2, but not at <f> = 0.0 as for the optical and ultraviolet peak 
(see §3.2). 



3.1. WFCAM Photometry During Minimum 



It is somewhat surprising that the nominal WFCAM H and K magnitu des in Table [2] 
are 0.12 — 0.15 mag fainter than the minimum reported by lDebes et al.l (120061 ). although the 
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datasets have overlapping error bars at each passband. As can be seen in Figure |3l it has 
been verified that the WFCAM observations were taken during the p hotometric min i mum b y 
comparing the corresponding corrected MJD with the ephemeris of iBurleigh et al.l (j2006al ). 



Five of the stars listed in Tabled] were used as calibrators by lDebes et al.l (120061 ). How- 
ever, the 2MASS photometric errors on 3 of these objects, while mainly decent at J, become 
relatively high at H and K s (erj < 11%, a H < 17%, ok s < 25%), and were discarded here 
for this reason. In order to investigate whether this fact could account for the discrepancy 
between the reported HK magnitudes - and more importantly, to assess the likelihood that 
this difference is real - the instrumental zeropoints for the 3 discarded calibrators were com- 
pared with those derived with the 12 higher S/N stars used here. Those 3 stars alone yield 
an H Q and K Q which differ by 3 standard deviations (—0.06 mag at H, —0.14 mag at K) 
from the mean zer opoints used for the Table [2] WFCAM results. Additionally, according to 



Debes et al.l (120061 ). not all 5 of these calibrators were present in every frame of their WIRC 



and PANIC data. The 3 stars with lowest S/N (those discarded here) are in fact closest 
on the sky to SDSS 1212, at 61 — 93", while the remaining 2 calibrators are 115 — 126" 
distant. Hence the 3 lower quality 2MASS calibrators had a good chance of being present in 
most, if not all, the WIRC frames (a single chip is 201" x 201"), while the 2 higher quality 
calibrator stars were unlikely to be imaged in the majority of that dataset. Owing to the 
small er field of view of P ANIC (128" x 128"), only the 3 lowest S/N 2MASS calibrators were 
used JDebes et alil2006h . 

There are two additional reasons why the WFCAM magnitudes in Table [2] may be 
slightly discrepant from the previously published values. First, neither study has transformed 
the photometry from the 2MASS system into the native filter system through which the 
observations were taken. In the present case, any errors introduced by failing to transform the 
photometry are much smaller than the errors due to the intrinsic S/N (represented by a p hot m 
Table [2]). Also, the fact that 12 calibrators were used here would a lmost certainly shrink any 
error of this type even further due to averaging (lCarpenterll200ll ). Second, it is conceivable 
that SDSS 1212 exhibits further variability which has so far not been documented. 



3.2. NIRI Photometry Near Maximum 



The NIRI photometry is also listed in Table [2J It has been verified that these observa- 
tions were taken near the minimum hot spot viewing angle. Although this corresponds to 
a maximum in ultraviolet and optical light, it does not quite correspond to the cyclotron 
maximum as seen in the K-band light curve measured by lDebes et al.l (120061 ). which occurs 
when the hot spot is on either limb - immediately before and after the ultraviolet maximum 
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0.0 in this paper and iBurleigh et al.ll2006al . <j) = 0-5 in lDebes et al.ll2006l ) 



3.3. GNIRS Phase- Averaged Spectroscopy 

The spectroscopy essentially covers all orbital phases of SDSS 1212 because the flux 
was sufficiently low (0.1 mJy at 1.6 /im) as to preclude extracting phase-resolved data - the 
observation itself represents a benchmark near the limit of what is possible with the most 
advanced ground-based telescopes and instruments. The average of all data from both nights 
yields S/N ratios in the range 10 — 17 (average of 14) at H and 10 — 26 (average of 18) at 
K, estimated using 12 sections of 50 pixels each over the spectral ranges plotted in Figure 
[2j Without any inter-order adjustment, the relative levels of the three shortest wavelength 
orders appear to match each other well in addition to being consistent with the J and H 
photometry. 

In order to perform a check on the flux level of the i^-band spectroscopy, the expected 
photometric flux in the K-benad was calculated based on the timing of the spectroscopic 
observations over the known X-band variability represented in Figure [31 Taking K m ^ n = 
17.37 mag from Table [2] and adding t o that the —0 . 70 m ag difference between the average 



.fT-band minimum and maximum from lDebes et al.l (120061 ). yields K max = 16.67 mag. From 



Figure [31 the spectroscopic observations span, over time, an estimated average minimum to 
maximum ratio of 5:7, yielding a flux-averaged magnitude of K avg = 17.02 mag over all the 
exposures. All three of these K-band fluxes are plotted in Figure [21 As can be seen in the 
figure, the spectroscopic flux over the X-band is commensurate with expectations without 
any adjustments. 

The spectrum in Figure [2] appears mostly consistent with a Rayleigh- Jeans slope in the 
J-band, where an emission line is likely detected at Paschen f3 and possibly also at Paschen 7. 
It is unclear whether any emission originates on the companion via heating by the primary, 
or on the white dwarf via accretion; both are possible. A modest excess is detected over 
the iZ-band, and is consistent with the Table [2] photometry. At ZT-band there is a definite 
excess throughout the window and four possible noteworthy features with varying degrees 
of uncertainty: a slight peak near 2.1 fxm that could be a cyclotron harmonic, Brackett 7 
emission, the CO bandhead at 2.29 /im, as well as a significant upturn in slope near 2.3 /im. 
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ANALYSIS 



4.1. Binary Components 

The near-infrared spectrum of SDSS 1212 has been fitted with several sp ectral models 



and te mplates. The binary components were analyzed following the method of iDobbie et al. 
( 120051 ). A model white dwarf spectrum spanning 0.3 — 10 /mi was genera ted for a non- 
magne tic, pure hydro gen atmosphere at T P ff = 10, 000 K and log g = 8.0 (ISchmidt et al. 
2005a ). using tlustyQ ( Hubeny Sz Lanz 1995 ) and synspec 1 . Since the white dwarf has a 



significant magnetic field, with Zeeman-split Balmer lines in its optical spectrum, it is not 
possible at present to more precisely determine its parameters by the usual method of Balmer 
line modeling. The model spectrum has been normalized to the g magnitude of the white 
dwarf and smoothed to match the spectral resolution of the GNIRS data. 

Template 1 — 2.3 /mi spectra for brown dwarfs of late L and early T type were taken from 
low resolution Keck / NIRSPEC datasets (IMcLean et al.ll2003l ). These templ ates were then 



extended to 2.5 /mi using similar near-infrared UKIRT / CGS4 spectral data (IGeballe et al. 



20021 ; iLeggett et al.ll200ll ). in order that they span the full extent of the GNIRS wavelength 
coverage. The NIRSPEC spectral templates have similar resolution (R « 2000) to the 
GNIRS data and hence no smoothing of either data was performed. However, the CGS4 
data had to be interpolated between resolution elements in order to create sufficient data 
points to compare with the GNIRS spe ctrum. Empirical model fluxes, determin ed from near- 
infrared photometry and parallaxes (jTinney. Burgasser. fc Kirkpatrickl 120031 ) were scaled 
appropriately to d ~ 150 pc, the nominal distance to a white dwarf with the assumed 
parameters listed above. 



4.2. Cyclotron Emission 



A physically realistic model for cyclo tron emission was calcu lated utilizing the formu- 
lation and methods described in detail by lHoard et al.l (120071 ) and iBrinkworth et al.l (120071 ). 
The model employs five critical parameters to produce the emergent cyclotron spectrum. 
The magnetic field strength and electron temperature determine the cyclotron harmonic 
wavelengths, while the viewing angle affects the overall shape of the peaks. Also included 
are a parameter related to the transition from optically thick to optically thin emission (typ- 
ically well-constrained by the observations), and a scaling factor linked to the area of the 
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emitting region; a few percent of the white dwarf surface area for the models shown here. 



Initial attempts to fi t the spectrum of SDSS 1212 using cyclotron emission in a B = 13 
MG ( ISchmidt et al.ll2003l ) field were problematic. The sum of the white dwarf and secondary 
model components does a decent job of reproducing the phase-averaged spectrum shortward 
of i^-band. Yet the addition of such a cyclotron model component in this bandpass does not 
account for the upturn near 2.3 /im in the GNIRS data. Several circumbinary dust models 
were considered to explain this observed rise in flux at the end of the JT-band, but were 
subsequently discarded upon the realization that this feature must be a cyclotron harmonic 
for the following reasons. 



As pointed out by lDebes et al.l (120061 ). the zero-temperature m = 4 cyclotron harmonic 
for a 13 MG magnetic field lies near 2.1 /im. But in such a model, this cyclotron feature 
alone, varied by 100%, is insufficient to produce the peak-to-peak flux change (Am « 1 mag) 
seen in their K-band light curve. Once the contributions of the white dwarf and secondary 
are accounted for, all the remaining flux in the .fT-band must be 100% variable in order to 
explain the observed light curve; leaving cyclotron as the only viable candidate for this extra 
emission. Furthermore, a 13 MG magnetic field should produce cyclotron harmonics and 
photometric variability wit hin the if -band (m = 5 would be near 1.6 /im), as seen in EF 
Eridani and AM Herculis (jKafka et al.l 120051 ; lHarrison et al.l 120041 ) . The resolution to this 
dilemma is that the magnetic field of SDSS 1212 must be less than 13 MG, such that more 
cyclotron emission is present within the K-band, while the harmonics at shorter wavelengths 
are diminished. 

A magnetic field of 7 MG should have two zero-temperature cyclotron ha rmonics within 
the K -band, near 2.2 and 2.5 /im for m = 6 and 7 respectively. Although ISchmidt et al. 
(120031 ) determine an equivalent dipolar field of B<i = 13 MG for SDSS 1212 (via models and 
assumed geometry), the mean surface field as determ ined by the measured Zeeman splitting 
of Ha and H/? is B s = 7 MG (ISchmidt et al.l l2005al ) . As shown below, this magnetic field 
strength is consistent with the absence of observable cyclotron harmonics in if-band, and the 
presence of two harmonics in the K-bernd which together can explain the large photometric 
variability seen in this filter. 

Figures [4] and [5] display the most successful model reproductions of the GNIRS cross- 
dispersed spectrum of SDSS 1212. The essential difference between the two figures is that 
an L8 dwarf template was used for the Figure H] plots, while a T2 dwarf template was used 
in the plots of Figure It should be noted that between these two figures, a minor scaling 
readjusment was made to the binary components in order to best reproduce the observed 
spectrum of SDSS 1212. This rescaling amounted to 8%, corresponding to a 4% change 
in the assumed distance to the system; well within the 14% minimum uncertainty arising 
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from its poorly constr ained temperature and unknown surface gravity (IBurleigh et al.ll2006a 



Schmidt et alJl2005ah . 



The upper panels in both figures represent the best fit to the spectrum with only the 
binary components at d ~ 150 pc; the white dwarf model plus a brown dwarf template. 
Already it is clear that the ff-band excess is largely reproduced by the addition of the L8 
dwarf, but not as well by the T2 dwarf. However, in both upper panels there appears to be 
excess flux over the 1.05 — 1.20 /iin range compared to the model. It is not clear if this is 
the result of poor signal or calibration in order 6 of the cross-dispersed GNIRS data, or a 
real detection. 

The middle panels of both figures add a single field cyclotron emission component to 
the upper panel model fits. In both figures, the models are for a 7.6 MG magnetic field, a 5 
keV plasma temperature, and a viewing angle of 75°, with the remaining parameters varied 
slightly so that the peaks of the cyclotron harmonics match the phase-averaged X-band flux. 
Although this set of models can account for a significant portion of the K-b&nd flux and 
hence go a long way towards accounting for the photometric variability in this window, they 
fail miserably in the region between the cyclotron peaks. It is perhaps noteworthy that the 
T2 model requires a stronger cyclotron component than does the L8 model (due to the bluer 
H — K color of T dwarfs versus L dwarfs) and is more readily consistent with the X-band 
variability. 

The lower panels in both figures display cyclotron models similar to the middle pan- 
els, but with a viewin g angle of 50°; a value commensurate with the determination by 



Burleigh et al.1 (l2006al ). Although this angle is not strongly constrained by their data, 
the GNIRS spectra appear to demand an intermediate viewing angle, whose result is to 
broaden the cyclot ron harmonics, effectively blurring them together for low viewing angles 



( iHoard et al.1 120071 ). In Figure HI the magnetic field employed is 7.0 MG , while in Figure 
[5] it is 7.2 MG, both adjusted slightly to best fit the overall spectrum when combined with 
their respective secondary models. Both these model fits do a decent job of reproducing the 
entire spectrum from 1.2 — 2.5 /xm, but the L8 model matches the H- and JT-band data more 
closely than the T2 model, as in the middle panels. 



4.3. Substellar Secondary 

The three component, composite model fits contain a sufficient number of free param- 
eters that their success in reproducing the 1.2 — 2.5 /im spectrum of SDSS 1212 may not 
seem surprising, but is nonetheless informative. During initial attempts to model the GNIRS 
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spectrum, secondary spectral templates of types L6 and T5 were shown to be significantly 
too bright and too dim, respectively, beyond 1.5 /im when their overall J-band flux level 
was matched to the data. Figure @] is consistent with the direct spectroscopic detection of a 
substellar object with spectral type L8. While the T2 model in Figure [5] appears somewhat 
inconsistent with the data, it does not rule out an earlier type T dwarf, for which available 
templates were lacking. Hence, the companion to SDSS 1212 has a likely spectral type be- 
tween L8 and Tl. Uncertainty in the distance to the white dwarf allows for some leeway 
in the relative contribution of the binary components at near-infrared wavelengths. While 
its poorly constrained effective temperature permits photometric distances both nearer or 
farther than d ss 150 pc, the unconstrained mass and radius of the white dwarf are likely to 
overpredict its distance. A typical isolated white dwarf surfac e gravity, or log g = 8.0 (0.60 
M Q at T eff = 10, 000 K; iBergeron. Saumon. fc Wesemaellll995l ). is more likely to be too low 
than too high, especially if the primary has been accreting material over Gyr timescales. 
This would tend to reduce the photometric distance and lead to later spectral types for the 
companion. If SDSS 1212 is more massive than log g = 8.0, there is a greater likelihood that 
the companion has an early T spectral type. The lower panel in the L8 model, which also 
matches the observed K-band spectrum well, seems to indicate that the CO bandhead at 
2.29 /im may have been detected, albeit just barely. If correct, it would be further evidence 
that the companion is a late L dwarf, since this CO feature weakens signific antly in T dwarfs 
JCushing. Ravner. fc Vaccall2005l : iMcLean et aDl2003l : beballe et al1l2002h . 



4.4. The Nature of SDSS 1212 



The two models for the nature of mass transfer and accretion in SDSS 1212 are Roche 
lobe overflow and wind capture. The near-infrared observations do not distinguish between 
these two evolutionary scenarios, but are informative for the following reason. The confirmed 
presence of cyclotron radiation adds one more observational property that SDSS 1212 now 
shares with polars (especially those with likely substellar secondaries in persistent low states 
such as EF Eridani). The traits it shares are: an appreciable x-ray luminosity, the presence 
of an accretion hot spot, a variable optical light curve with a similar shape to those of known 
polars in low states (i.e. with larger modulations in the blue and ul traviolet compared to 
the red), a spin- locke d binary configuration, and cyclotron emission (IBurleigh et al.ll2006al ; 
Schmidt et al.ll2005al ). Occam's razor would argue that SDSS 1212 is a polar. 

Furthermore, the hot spot x-ray luminosity of SDSS 1212 is L x « 3 x 10 29 ergs s , which 
is nearly i dentical to that meas ured for EF Eridani during its 10 yr low; namely L x « 2 x 10 29 



ergs s 1 (jSchwope et al.l 120071 ). While the x-ray emission from many suspected, detached, 
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wind capture systems may be too weak (and more importantly, not periodically variable) 
to originate via accretion onto the white dwarf, but rather from activity on the M dwarf 
secondary, L dwarfs simply do no t have sufficient x-ray luminosities, hence favoring the low 
state polar model for SDSS 1212 (IBurleigh et al.ll2006al ). Also, there is the issue of how such 
a low mass object generates sufficient wind to be captured by the magnetic white dwarf. To 
date, no model has emerged which can explain how an old L or T dwarf might generate a 
wind on the order of 10~ 13 M Q yr _1 . 



5. DISCUSSION 
5.1. Past and Future Evolution 

The gross evolution cataclysmic variables should consist of three fundamental phases: 

(1) A short- to long-lived phase as a detached main sequence binary which subsequently 
evolves to shorter orbital periods through common envelope evolution as the primary leaves 
the main sequence and becomes a first ascent or asymptotic giant. Angu lar momentum i s 



transferred, via friction, into the envelope itself, ejecting it in the process (jPaczynskil Il976l ) . 
If the orbit is still too large for the secondary to make contact with its Roche lobe, then three 
mechanisms can bring this about, in principle: gravitational radiation, magnetic braking, or 
nuclear evolution. The timescales of these processes depend on the masses involved, but if 
one assumes a canonical 0.6 M & white dwarf and a 0.2 M secondary (corresponding to the 
mass distribution peak among detac hed, unevolved, low mass companions to white dwarfs; 



Farihi. Becklin. fc Zuckermanll2005al ). then both gravitational radiation and magnetic brak- 



ing operate within Gyr timescales for period ranges P < 3 hr a nd P > 3 hr res pectively. 



whereas nuclear evolution will play no role within this timeframe ([Patterson! 119841 ) . 



(2) A phase in which the shrinking Roche lobe of the secondary becomes smaller than its 
radius, causing unstable mass transfer onto the white dwarf. In principle this is a runaway 
inspiral until the mass ratio is close to one, but if the secondary has M < 0.5 M & to begin 
with (as it should in the vast majority of cases based on the aforementioned companion mass 
distribution), then conservation of angular momentum competes to partiall y counterbalance 



the lo sses due to any braking (magnetic or tidal) and gravitational radiation (INelemans et al. 



200ll ). At some point the trend reverses outward as the orbital evolution becomes dominated 
by conservation of angular momentum and any mass lost from the system entirely (e.g. nova- 
like outbursts, accretion-driven winds). 

(3) A final, debatable, phase consisting of a once more detached or very low accretion 
rate semi-detached binary comprised of two degenerate objects in an expanded orbit which 
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persists (now subject only to gravitational radiation) up to 10 Gyr or so, depending on 
the final masses and separation. This eventual state depends on the thermal and nuclear 
history of the secondary but requires a degenerate secondary core in order to cease expansion 
permanently. Thus, a system may end its existence as a white dwarf with a close, altered 
substellar companion, both objects cooling gradually toward oblivion. 

The ultimate fate of catacylsmic variables outlined above is uncertain for two funda- 
mentally different reasons. The first issue is that the outcome itself is unknown, with end 
predictions for secondaries including; complete stripping, tidal disruption, persistence as 



substellar or even planetary mass objects ( 


Littlefair. D 


lillon. & Martin 


2003; 


Kolb & Baraffe 


1999; 


Patterson 


1998; 


Howell. Rappaport. & Politano 


1997 


). The second issue is detectabil- 



ity, which itself divides into two categories: searching for either erupting systems or their 
burnt out remnants. 



5.2. Recently Deceased or Born Dead? 



So far, astronomers who study cataclysmic variables have concentrated primarily on 
those systems known to be eruptive in order to search for and understand possible end states, 
including substellar companions. The means of identifying such low mass secondary systems 
relies heavily on semi-empirical and theoretical relations between obser ved periods and spec 



tral type, mass ratio, effective temperature, and mass transfer rate (ILittlefair et al.l 2006 
20071: Harrison et"aUl20od:IPatterson. Thorstensen. fc KemrJ[2005l: iKolb fc Baraffe 1999 



Howell. Rappapor 



19971 ). Without models, there is no means of determining whether a given substellar com- 
panion is migrating inward or outward, i.e. approaching or leaving its period minimum. 
With the prior history of these interacting systems unknown, the beginning and end states 
may appear identical, theoretically. 

Additionally, it is difficult to know with certainty what the size of the seconday star 
is in such systems, especially if the mass transfer rate is low - consistent with quiescence 
or an underfilled Roche lobe (i.e. a non-catacylsmic variable). If a companion does indeed 
have a radius smaller than its Roche lobe, then strictly speaking it is detached, and the 
question remains: is Roche lobe contact in its f uture or in its past? T he recently identified 
low mass transfer magnetic binaries reported in ISchmidt et al.l (j2005bl ) are good candidates 
for being detached due to the combination of the inferred accretion rates, orbital periods, 
and extremely low or absent x- ray emission. But if EF Eridani can turn off for nearly 
10 yr, then so can other polars (jSchwope et al.l 120071 ) . SDSS 1212 sits somewhere on the 
border between a typical polar and the low accretion rate systems, having respectable x-ray 
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accretion luminosity and a relatively low inferred mass transfer rate (IBurleigh et al.ll2006al ). 
If one assumes the substellar companion to SP SS 1212 has a mass an d effective temperature 
similar to WD 0137-349B (L8V, see Table 131 liurleigh et al.ll2006bh then its unperturbed 
radius should be 0.09 R Q for ages between 1 — 5 Gyr (IBaraffe et al.l 120031 ) . but its Roche 
lobe should be located roughly at 0.11 Rq for a 0.60 M & primary. Therefore, the companion 
is either: 1) recently deceased, far out of thermal equilibrium and inflated as expected for 
an object past period minimum with the implied long history of mass-radius imbalance; or 
2) born dead, intrinsically detached, and slamming the magnetic white dwarf with all of its 
(irradiatively driven?) wind. This raises similar questions about the nature of all other low 
accretion rate magnetic binaries. 

Depending on the past thermal and nuclear history of the secondary, especially if the 
core is degenerate, and how much mass escapes the binary completely and adiabatically, 
it is quite possible for the orbital period to increase signifi cantly over the minimum into a 
region where no i nteraction (or wind accretion only) occurs ( jHowell. Rappaport. fc Politano 
19971 ; |Jeanslll924l ). Without any distinct spectral features or mass-radius anomalies in such 
detached binaries, there may be little or no means to differentiate empirically between the 
pre- and post- states of these innately interacting systems. 



5.3. The Frequency of the Doubly Dead 



Objects must first be detected as a catacylsmic variable in order to look for substellar 
companions; hence, selection effects play a l arge role in the number of system s suspected 
to harbor such cool degenerate companions (ILittlefair. Dhillon. fc Martini 120031 ) . However, 
this overlooks the obvious question: if the end states are detached double degenerate sys- 
tems, where are they? If one restricts a search to eruptive systems, then the potential for 
success is limited from the beginning for reasons stated above. Why not search directly for 
systems which are detached and which may represent final states? This seems to be the 
only way to reconcile the well-known discrepany between the number of p redicted objects 



past period bounce and the number o bserved or suspected from observation (jPattersonlll998 
Howell. Rappaport. fc Politano! 1 19971 ). 



Clearly, the best place to look for these types of end states is in the i nfrared where cool, 

Jupiter-sized companions easily outsh i ne thei r tiny white dwarf primaries (IFarihi. Becklin. fc Zuckerman 
2005al : IFarihi. Zuckerman. fc Becklin! l2005bl ) . There are only two issues associated with a 
search around white dwarfs for extinct cataclysmics - the white dwarf must first have been 
identified or previously known and one must have sufficient sensitivity to detect its cooling 
leftovers. If the end state white dwarfs are very cool and far away such that V > 21 mag 
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(e.g. T eff < 5000 K and d > 100 pc for M = 0.9 M ), then only the Sloan Digital Sky Survey 
will detect them but will be quite insensitive to any substellar companions. While the above 
conditions may not be representative, given that cataclysm ic variables appear to be roughly 
100 times less abundant than white dwarfs (jPattersonl 1 1 984l ) . a large sample of targets would 
be required to perform a statistically rigorous search; if all cataclysmic variables are in this 
state, a search would only yield 1 per 100 white dwarfs (not very encouraging). 

If these close but detached systems exist in any significant numbers, they have not yet 
been identified from various surveys of relatively nearby white dwarfs. Spatially unresolved 
L type companions, corresponding to 0.05 — 0.07 M & for unevolved ages of 1 — 5 Gyr, are 



rare, with only 3 known among well ov er N = 1000 1 targets Jlloard et al.ll2007l: 



Burleigh et al. 



2006bl ; iFarihi. Becklin. fc Zuckermanll2005al : iFarihilhooi iFarihi fc Christophe rlbooi IvVachter et al 
20031 ). Spatially unresolved T type companions, corresponding to 0.02 — 0.05 M for un- 
evolved ages of 1 — 5 Gyr, are also rare, with none detected or suspected in various sur- 
veys totalling N 200 targets (jMullalv et al.ll2007l ; Hansen. Kulkarni. fc Wiktorowiczll2006l ; 



Farihi. Becklin. fc Zuckerman 2005al ; Farihi. Becklin. & Zuckerman 2008, in preparation). 
It is possible that stellar evolution has yet to produce such systems; i.e. the initial mass 
function for secondaries in such systems, convolved with the finite age of the Galaxy, implies 
that the vast majority of these systems will continue in phases (1) and (2) for up to 10 Gyr. 
It could be that the lack of detected end states is an observational bias; if these systems 
are still marginally semi-detached then the signature of the secondary star could be masked 
by the accretion disk in dwarf novae and by cyclotron emission in polars. Alternatively, 
phase (3) never occurs and the companion is cannibalized, evaporated, stripped or tidally 
destroyed. This would also explain the lack of observed end states. 

For completeness, something ought to be said ab out the empiric al frequency of zero age 
cataclysmic variables with a brown dwarf secondary (lPolitanoll2004l ). Since their pre- states 
would consist of a white dwarf with a detached substellar companion, the above arguments 
and evidence apply (/ < 0.5%). Additionally, the progenitors of such pre- states would 
be main sequence stars with brown dwarf companions in the few AU ra nge, which are also 



extremely rare, oc curing with a well-constrained frequency, / < 0.5% (IButler et al.l 12006 



Marcy et al.ll2005l ). Hence, there should be no doubt that cataclysmic variables form with 
brown dwarf secondaries quite rarely. 



5.4. Detached Companions to Magnetic White Dwarfs 



The same argument which applies to detached substellar companions to white dwarfs 
also applies to detached stellar companions to magnetic white dwarfs, although in this latter 
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case finding them should be easier in principle (ILiebert et al.ll2005l ). One might suggest a 
sensitive search for magnetism among known lists of close but detached, and wide common 
proper motion companions. Large lists of these white dwar f plus red dwarf systems exi st in 
the following references: for known radial veloc it y pairs see 
spatially unresolved pairs see Hoard et al.l (120071 ) ; iFarihi. Hoard, fc Wachter 
(2005a); 1- 
j2005ah : 



Wachter et al.l (120031). and for wide common proper motion pairs see 



d2005j), for 



2006 N ) :lFarihi. Becklin. fc Zuc s 



Farihi. Becklin. fc Zuckerman 



Silvestri. Hawley. fc Oswaltl (120051 ). All of these lists contain targets which are likely 



to be superior to white dwarfs from the Sloan Digital Sky survey since they are brighter and, 
hence, more sensitive searches could be conducted. However, it should be noted that the 
main observational indicators of a magnetic white dwarf, Zeeman splitting and circular spec- 
tropolarimetry, are strongest in Ha and H(3. Unfortunately, in unresolved white dwarf plus 
red dwarf binaries (i.e., the potential cataclysmic variables) these features, especially Ha, 
are often completely filled-in by the red dwarf. This may be another reason for the curious 
lack of detected pre -polar systems to be considered along with those already discussed by 



Liebert et all (120051 ). 



6. FUTURE WORK 

Unfortunately, the intrinsic faintness of SDSS 1212 precludes several interesting avenues 
of follow up, such as phase-resolved near-infrared spectroscopy or mid-infrared spectroscopy. 
Ultraviolet and x-ray observations are currently planned to firmly distinguish between stan- 
dard mass transfer through Roche lobe overflow versus accreted wind. High resolution optical 
spectra of SDSS 1212 may reveal the presence of elements accreted from its companion. Such 
observations could, in principle, distinguish between Roche lobe overflow mass transfer and 
wind accretion as well. 

There may be little or no hope of directly detecting the bare substellar companion be- 
yond the i?-band spectrum presented here. The 7 MG cyclotron models used successfully 
in this paper predict the fundamental harmonic wavelength should be 15.3 /mi, thus po- 
tentially swamping every instrument aboard Spitzer which might otherwise be capable of 
directly detecting the flux of an L8 dwarf at d = 150 pc; namely IRAC 3 — 8 /mi photome- 
try, and IRS 16 /mi peak-up imaging. However, IRAC observations should directly test the 
lower magnetic field cyclotron models employed to fit the GNIRS data, which predict rela- 
tively strong emission near 3.1, 3.8, 5.1, and 7.7 /mi, corresponding to the zero-temperature 
m — 5, 4, 3, and 2 harmonic numbers in a 7.0 MG magnetic field. Yet the poss ible presence 



of circumbinary dust emitting at IRAC waveleng ths, as seen in several polars (IHoard et al. 



20071 ; iBrinkworth et al.l 120071 : iHowell et al.l 120061 ) . may complicate the interpretation of any 
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mid-infrared photometry for SDSS 1212. 
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Table 1. WFCAM Photometric Calibration Data 



Number 


2 MASS Name 




H (mag 


) 


K s (mag 


) 


#o f ( 


mag) 


( 


mag) 


Remarks 


1 


J12114360+0137020 


13 


841 


±0 


035 


13.701 ±0 


050 


23 


032 


22 


354 




2 


J12114583+0139374 


13 


435 


±0 


030 


13.337 ±0 


038 


23 


083 


22 


391 




3 


J12114736+0133594 


14 


098 


±0 


039 


14.105 ±0 


063 


23 


031 


22 


426 




4 


J12114758+0137223 


14 


021 


±0 


037 


14.062 ± 


062 


23 


040 


22 


479 




5 


J12115403+0132049 


13 


991 


±0 


046 


13.766 ± 


058 


23 


069 


22 


325 




6 


J12115541+0130399 


13 


431 


± 


027 


13.299 ±0 


035 


23 


041 


22 


356 




7 


J12120408+0135365 


15 


046 


±0 


066 


14.776 ± 


123 


23 


076 


22 


282 


1,2 


8 


J12120582+0135157 


14 


987 


±0 


072 


14.955 ± 


138 


22 


980 


22 


404 


1,2 


9 


J12120662+0138153 


14 


065 


±0 


053 


13.974 ±0 


068 


23 


075 


22 


405 




10 


J12120989+0135259 


15 


894 


±0 


174 


15.538 ±0 


232 


22 


908 


22 


115 


1,2 


11 


J12121658+0137042 


14 


290 


±0 


041 


14.137 ±0 


066 


23 


039 


22 


452 


2 


12 


J12121667+0135257 


13 


868 


±0 


022 


13.805 ±0 


068 


23 


036 


22 


383 


2 


13 


J12122238+0130101 


12 


927 


±0 


022 


12.904 ±0 


032 


23 


058 


22 


410 




14 


J12122266+0137425 


13 


602 


±0 


041 


13.536 ±0 


055 


23 


Oil 


22 


470 




15 


J12122647+0140477 


13 


040 


±0 


029 


13.030 ±0 


037 


23 


070 


22 


412 





Note. — (1) Not used in the present paper; (2) Used in iDebes et al.l (120061 ) 
''"Transformations between 2MASS and UKIRT photometric systems were ignored (see §3.1). 
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Table 2. Photometry of SDSS 1212 



Passband 


Instrument 


Magnitude (mag) 


a pho t (mag) 


o ca i (mag) 


atotai (mag) 


H 


WFCAM 


17.676 


0.089 


0.021 


0.091 


H 


NIRI 


17.710 


0.064 


0.030 


0.071 


K 


WFCAM 


17.372 


0.105 


0.047 


0.114 


K 


NIRI 


16.693 


0.033 


0.030 


0.045 



Note. - An intrinsic calibration uncertainty of 3% was assumed for the NIRI data. 
See Figure [1] for the phases corresponding to the photometric observations. 
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Table 3. Close Binary Parameters for SDSS 1212 and WD 0137-349 



Star 


P (min) 


M 1 (M Q ) 


M 2 (M ) 


a (R & ) 


R 2 (Rq) 


Rl2 (Re) 


WD 0137-349 


115.6 


0.39 


0.05 


0.60 


0.09 


0.13 


SDSS 1212 


88.4 


0.60 


0.05 


0.57 


0.09 


0.11 



Note. - Values for the white dwarf and secondary masses in the SDSS 1212 sys- 
tem are assumed since no measurements exist. However, the resulting semimajor axis, 
unperturbed secondary radius, and secondary Roche lobe radius are all relatively in- 
sensitive to those input values, and reflect likely parameters. The semimajor axis for 
WD 0137—349 was ca lculated using Kepler 's law, yielding a value which differs slightly 
from that reported in iMaxted et al.l (120061 ). 
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Fig. 1. — Spectral energy distribution of SDSS 1212. Optical gri and near-infrared zJ 
photometry is from ISchmidt et al.l (l2005al ). The i?-band value is the average of the NIRI 
and WFCAM measurements, while the K-band value is the flux measured during minimum 
(WFCAM), all from the present paper (Table [2]). The grizJ photometric fluxes are well 
reproduced by a 10,500 K blackbody. 
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Fig. 2. — Flux calibrated GNIRS spectrum of SDSS 1212, averaged over all nights and 
phases. Four cross-dispersed orders are displayed within the atmospheric windows. The 
photometric points and blackbody continuum are identical to that plotted in Figur^H with 
2 additional (calculated) i^-band flux points described in §3.3. 
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NIRI 
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Fig. 3. — Orbital phases corresponding to the photometric and spectroscopic observations 
with each instrument: GNIRS, NIRI, and WFCAM. The numbers 1 — 6 denote each of 6 
nodded pairs of GNIRS observations taken on two separate nights (subscr ipted 1 and 2). 
Also shown are the it-b and light curve m odel fit of iBurleigh et al.l (j2006al ) and the rough 
.fT-band light behavior flDebes et al.l 120061 ). both normalized to one. The relative amplitude 
of the ti-band peak has been exaggerated by a factor of 10 for visibility. 
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Fig. 4. — Model fits to the GNIRS, phase-averaged, near-infrared spectrum of SDSS 1212. In 
each panel the data are represented by a dotted line and the fit is shown as a solid line. Some 
rescaling has taken place between panels to best fit the data, but the binary components 
remain at d ph 150 pc in all plots. Details of the models and fits are discussed in §4.1 and 
4.2. Top panel: The dashed line is a pure hydrogen atmosphere, T e g = 10, 000 K model for 
a 0.6 M & white dwarf. The dash dot line is an L8 brown dwarf template spectrum. Middle 
Panel: The dash dot dot line shows the cyclotron model at a viewing angle of 75°. Bottom 
panel: Here, the dash dot dot line shows the cyclotron model at a viewing angle of 50°. 
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